Longitudinal spin-fluct nations and superconductivity in ferromagnetic ZrZn2 from 

ab initio calculations 
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The recent discovery of superconductivity coexisting with 
weak itinerant ferromagnetism in the d-electron intermetal- 
lic compound ZrZn2 strongly suggests spin-fiuctuation me- 
diated superconductivity. Ab initio electronic structure cal- 
culations of the Fermi surface and generalized susceptibilities 
are performed to investigate the viability of longitudinal spin- 
fluctuation-induced spin-triplet superconductivity in the fer- 
romagnetic state. The critical temperature is estimated to be 
of the order of 1 K. Additionally, it is shown that in spite of a 
strong electron-phonon coupling (Aph = 0.7), conventional s- 
wave superconductivity is inhibited by the presence of strong 
spin-fluctuations. 



The generalization of the Bardeen-Cooper-SchriefFer 
(BCS) theory to electron-electron interactions by Kohn 
and Luttinger |l|, paved the way for speculation about 
the possibility of non-s-wave, or "unconventional" super- 
conductivity. Following the suggestion that a magneti- 
cally mediated interaction plays an important role in the 
superfluidity of liquid "^He [p|, the search began for su- 
perconductivity in nearly magnetic metals where strong 
spin fluctuations might provide the pairing mechanism 
p,E|. Recent experiments on Sr2Ru04 p] have made it 
a strong candidate for exhibiting spin-triplet, possibly p- 
wave superconductivity. For a spin-singlet Cooper pair, 
where the electrons have anti-parallel spins, the presence 
of ferromagnetic spin-fluctuations will be antagonistic to- 
wards the development of such a superconducting state. 
However, the recent reports of the coexistence of ferro- 
magnetism with superconductivity in UGe2 [p| and ZrZn2 
[[7| suggest a spin-triplet Cooper pairing, probably driven 
by such spin-fluctuations. Moreover, in ZrZn2, the dis- 
appearance of superconductivity at the same point as 
magnetism, and the sensitivity of its occurrence to sam- 
ple purity [[7|,^ are perhaps the strongest indications yet 
that the superconductivity is intimately connected with 
the magnetism in this material. 

Unlike other "magnetic" superconductors (e.g. boro- 
carbides ||^ , RuSr2GdCu2 08 [1^) where the magnetism 
and superconductivity occur in different parts of the unit 
cell, in both UGe2 and ZrZn2 it is the same itinerant elec- 
trons that arc thought to form the Cooper pairs as well as 
produce ferromagnetism. Moreover, whereas some ques- 



tions regarding the itineracy of 5/ electrons and the roles 
of the strong magnetocrystalline anisotropy and quasi- 
two-dimensional electronic structure can be raised with 
respect to UGe2 , ZrZn2 is a three-dimensional intermetal- 
lic compound free of such effects. Discovered by Matthias 
and Bozorth ||l^ in the 1950s, it was initially of interest 
because of the presence of weak ferromagnetism, in spite 
of the fact that neither constituent was ferromagnetic. 
It has the C15 cubic Laves crystal structure, with a lat- 
tice constant of 7.393A (13.97 a.u.) [§. The possibility of 
there being triplet pairing in high-purity C15 compounds 
like TiBe2 and ZrZn2 was first suggested by Fay and Ap- 
pel [Q . In this Letter we investigate the viability of their 
suggestion in the case of ZrZn2. 

Firstly, we have calculated the electronic structure 
of ZrZna using the LMTO method |l|]. Exchange- 
correlation effects are described within the local spin 
density approximation (LSDA). Self-consistency was at- 
tained using 505 k-points within the irreducible wedge 
of the face-centered cubic Brillouin zone (BZ). The ba- 
sis included s, p, d and / states for all atoms. Our re- 
sults agree with the calculations of Jarlborg and Freeman 
]l3[ and Huang et al. p^ i n the non-magnetic state and 
those of de Groot et al. ]15[ and Jarlborg et al. M& in the 
spin-polarized one. The Fermi surface (FS), comprising 
four sheets, is shown in Fig. H. Both non-magnetic (NM) 
and ferromagnetic (FM) calculations were performed at 
a series of different lattice parameters, and the magnetic 
moment was found to disappear near 13.47 a.u., where 
the calculated pressure is 45 kbar. This is just below 
the total energy minimum, indicating a calculated equi- 
librium lattice constant of about 13.6 a.u., in excellent 
agreement with the recent FLAPW calculation of Bruno 
et al. [^. This underestimation (^2.5%) is typical of the 
LDA, particularly when including / states. Experiments 
confirm that the system is near the ferromagnetic insta- 
bility, since the critical pressure for the disappearance of 
magnetism has been reported in the range 8.5 — 22 kbar 

Polarized neutron studies of the magnetization density 
]l9[ have shown that there is a significant spin-density 
along the Zr-Zr bond directions. Our calculation also 
shows this delocalization, but typically with 85% of the 
total moment on the Zr. The dominance of Zr is also re- 
flected in the DOS at Ep, of which 70% is of Zr-d charac- 
ter. Mattocks and Dixon eO] inferred an exchange split- 



ting of 4.5 mRy (in a field of 8T) from their de Haas-van 
Alphen data for orbits on the F-centered spheroid (band 
30). This compares favorably with the value obtained 
from the current calculation at the equilibrium lattice 
constant (see Table |). 

Secondly, we have studied the FS and its nesting prop- 
erties. As seen in Fig. nl it is rich in details. Four bands 
(27-30) cross the Fermi level. We used the relaxed lattice 
parameter given by LSDA (a = 13.6 a.u.) corresponding 
to a total moment of 0.17 /is/Zr. All these sheets show 
a strong Zr-d character, though bands 27 and 30 exhibit 
a significant hybridization with Zn-p (25% and 30% re- 
spectively). Furthermore, the DOS at Ep is dominated 
by the contribution from band 29 (50%) and band 28 
(32%), whereas band 30 contributes less than 1%. The ef- 
fect of the spin-splitting is most noticeable in the change 
in the topology of the band 29 sheet, where the neck near 
the L point, present for majority spin electrons (f), dis- 
appears for minority ones (|). From Fig. |l|, we expect 
strong intra and inter-band nesting features, especially 
along the <100> direction. In order to understand how 
the nesting will actually affect the response of the elec- 
trons in this system, we have calculated the bare-band 
generalized static susceptibility. 
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where n denotes the band index, a the spin, and fank 
the Fermi-Dirac functions. This sum was calculated on 
a mesh of more than 5 x 10^ k-points in the cube (shown 
in Fig. n^) using a tetrahedron interpolation technique 
similar to that of Rath and Freeman [£1| . These calcula- 
tions were performed for both NM and FM cases along 
<100>, <110> and<lll>inq-space, at the LSDA equi- 
librium lattice constant. All these are peaked a,t q — 
and (rather surprisingly, given the FS topology) show 
very little structure a.t q > 0. This was expected for the 
NM case because of the ferromagnetism of ZrZn2, but the 
absence of finite q peaks in the FM case shows that this 
compound docs not favor antiferromagnetic (AF) spin- 
waves as confirmed by our frozen spinwave calculations. 
In fact, the only significant q > peak is found in the 
<100> direction at q = 0.2 (in units of 27r/a) in Xo (l)- 
This peak, originating from intra-band contributions in 
band 29 near the corner of the cuboidal FS sheet (i.e. 
near the L point) is (obviously) not present in the NM 
case. The importance of this peak (and correspondingly 
the absence of the other expected ones) is due to the con- 
centration of the DOS at Ep near the border of the BZ 
(points X, K and L) where band 29 flattens considerably. 
In other words, this means that while nesting is present 
elsewhere, it is inhibited by the low DOS. 

We now turn to the question of longitudinal spin- 
fluctuation-driven superconductivity in ZrZn2 as pro- 
posed by Fay and Appel H. On the FM side of the 



transition, when the band structure is different for the 
two spins, we calculate the longitudinal coupling con- 
stant, A^£, from the generalized susceptibilities. In the 
FM region, within the RPA, the pairing potential can be 
written as H 
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In the approximation of a spherical FS (not unreason- 
able, given the dominant influence of band 29 and its FS 
topology), the longitudinal coupling parameter is given 

byi, 
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where Pi is the Legendre polynomial. The objective is 
now to estimate the s and p components of Agf from 
our band structure calculation. The exchange integral, 
/, is obtained through the calculation of the Stoner en- 
hancement, S, defined as the increase of the exchange 
splitting of the Zr potential divided by the energy of 
the applied magnetic field. The corresponding Stoner 
factor 5^ = 1 — 1/ S is simply related to / through 
/ = S/N, where N is the DOS at Ep. We have cal- 
culated S{q) for the FM case (see Table Q) and for two 
AF spinwaves {q = Q and q = 27r/a). We find that S 
is quickly suppressed {S{q = 2-k / a) < 1.5) for AF spin- 
waves. This shows that ZrZn2 does not support AF fluc- 
tuations (which is consistent with the absence of peaks 
for g > in our calculated generalized susceptibilities) 
and that the ^-dependence of / cannot be neglected. 
From S{q), we model I{q) = /o/(l + b^q^) |,§ with 
/o = 0.04 Ry and b'^ = 0.33 (a/27r)2. Since thc_contribu- 
tion to Agf from each spinwave mode is ^S{q)S{q)'^ [E3|, 
we get, as our first estimate, Agf = 1.2 by averaging over 
these three modes. We make the further approximation 
Xo"^!?) = Ncr{Ep) Vg, and calculate the longitudinal X]^fi 
in the s {I = 0) and p {I = 1) channels from Eqs. (0) and 
(0), taking a kp„ « 0.6 27r/a, appropriate to the band 29 
sheet. For the relaxed FM case in which S — 1.12, A^j is 
negligible in the p channel. However, moving closer to the 
FM transition, i.e. for S = 1.01, we get A^/^p = 1.9,2.0 
(for a =t,i) in the s channel and A^/^^ = 0.81,0.76 (for 
a =t,i) in the p channel. This shift is justifiable given 
that S is extremely sensitive close to the transition. As 
noted by Fay and Appel, the s component is much larger 
than the p one and both diverge when S —^ 1. These 
values are consistent with our previous estimate. 

The electron-phonon interaction cannot be ignored and 
can even be expected to be rather large owing to the 
occurrence of conventional superconductivity in both Zr 
and Zn, and the large DOS at Ep. This suggests that 
the electron-phonon coupling could be sufficient to over- 
come the pair-breaking effects due to spin-fluctuations. 



The electron-phonon coupling constant, Aph, can be ex- 
pressed as Aph = X^i mTZJ^' where the sum runs over 
all atoms, i, with masses. Mi, and phonon frequencies, 
uji, while the numerator, ?/; = Ni{Ep){{WVi)'^), is the 
Hopfield parameter that describes the electronic contri- 
bution p4| . Here, 77 was calculated in the rigid muffin-tin 
approximation p4| , i.e. retaining only dipolar terms and 
neglecting electronic screening of the ionic displacements. 
The values for (cuf) are taken as one half of the Debye fre- 
quency of the atom i. Furthermore, we assume that the 
volume dependence of (ujf) follows VaB, where B is the 
(calculated) bulk modulus and a the lattice parameter, 
which is reasonable as long as all small-g phonons behave 
identically with pressure [E5|. As shown in Table |, Aph 
is of the order 0.7 near the calculated equilibrium volume 
and twice as large at the experimental lattice constant. 
Ignoring completely the destructive effects of spin fluc- 
tuations, and using /i* = 0.13 in the McMillan formula 
p6| , these correspond to respectable Tc's of about 8 and 
21 K, at the respective lattice constants. The decrease 
of Aph in the FM state can be attributed to the smaller 
total DOS at Ep and its pressure dependence can be as- 
cribed almost entirely to the behavior of the Debye fre- 
quency. In the vicinity of the ferromagnetic transition, 
this Aph is nevertheless insufficient to overcome the domi- 
nance of the spin-fluctuations as indicated by the large S 
enhancements near the critical pressure (see Table ||). A 
possibility exists for phonon-mediated superconductivity 
at larger pressures, i.e. well outside of the magnetic re- 
gion when the Stoner factor would be further decreased 
to completely suppress spin fluctuations. However, our 
calculations at a = 13.17 a.u. (equivalent to 160 kbar) 
indicate that Aph drops to 0.4 while S is still 4.1 (giv- 
ing a Asf of the order of 0.4), which combine to make 
the conditions non-favorable for phonon-mediated super- 
conductivity. Note that the persistence of such a large 
Stoner enhancement over this range of pressures shows 
again the importance of the spin fluctuations and that the 
large values for S in the FM region can be related to the 
observed absence of saturation of the magnetic moment 
[0 and the weak ferromagnetism. Since the magnetic 
moment, Stoner enhancement and spin fluctuations are 
associated with the Zr sublattice, phonon-mediated su- 
perconductivity might be envisaged to take place within 
the Zn sublattice, but such an explanation can be ruled 
out because the Zn contribution to rji is negligible. 

Having presented evidence against the possibility of 
electron-phonon driven superconductivity, we now try to 
estimate Tc from the longitudinal spin fluctuations. The 
typical spin-fluctuation cut-off frequency, ojgf , can be es- 
timated from the Stoner factor by Wgf — l/i^NS) p2[ , 
giving about 90 K at the relaxed lattice parameter. Us- 
ing the Allen-Dynes formula, we arrive at a simplified 
expression for the superconducting transition tempera- 
ture : 
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Note that the rather strong electron-phonon interac- 
tion Aph contributes to the mass renormalization (nu- 
merator) and is detrimental to superconductivity in this 
case. Furthermore, the s-wave A^f q contains both the 
longitudinal (L) and transverse (T) contributions. From 
the large measured electronic specific heat coefficient 
7cxp — 47 mJ mol^^K^^ 0], and our calculated values 
for A^/^Q and Aph, we infer a transverse contribution A J q 
of about 0.8. Using the values for the case S — 1.01, 
we get a T^ ~ 1.0, 0.8 K for a =|, I respectively. These 
estimates are very approximate but they confirm, in our 
opinion, the viability of triplet p-wave superconductivity 
in ZrZn2. 

In conclusion, we have shown that calculations based 
on our electronic structure results strongly support the 
idea that the recently observed superconductivity in 
ZrZn2 M is indeed a result of triplet pairing, as sug- 
gested by Fay and Appel 0] . This would lead to p-wave 
superconductivity, and since impurity scattering acts as 
a pair-breaker for pairing in the I 7^ channels, the high 
purity of samples is crucial. However, the experimental 
absence of superconductivity in the paramagnetic phase 
just above the critical pressure M is still unanswered by 
this theory which predicts an even larger Tc in the NM 
region. The answer may lie in the peak of the transverse 
susceptibility, Xo (9ioo), at q ~ 0.2 which could provide 
an attractive coupling that would naturally disappear 
outside of the FM phase. Finally, it might be worth- 
while revisiting the properties of C15 compound TiBe2 
under pressure since its electronic structure is very simi- 
lar to that of ZrZn2 and conventional superconductivity 
would be more favored owing to the lighter masses of its 
constituents. 
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TABLE I. Calculated parameters (per formula unit) for 
various lattice constants, a. Shown are the magnetic moment 
/i, the exchange splitting, ^, the Stoner factor, S, the density 
of states at the Fermi level for the non-magnetic (NM) calcu- 
lations as well as for the ferromagnetic (FM) ones in paren- 
theses (T/i), Debye temperature (set to 370 K at a = 13.573 
a.u.) used to calculate the electron-phonon coupling, Aph, in 
both the NM and FM cases, and the specific heat coefficient 
renormalized by the electron-phonon interaction, 7. 
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K (NM) (FM) mJ/molK" 



13.970 0.48 19.5 2.9 68 (18/27) 265 1.42 0.90 28.6 

13.573 0.10 5.0 8.3 54 (31/25) 370 0.71 0.72 16.1 

13.437 0.00 0.0 9.0 52 (26/26) 420 0.56 - 13.8 




FIG. 1. The spin-polarized Fermi surface for a = 13.6 
a.u., from bands 27 — 30 (top to bottom). The majority spin 
(t) sheets are shown on the left hand side and the minority 
spin (X) on the right. 






